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Abstract 

The PeV gamma ray background produced in the interactions of 
ultra high energy cosmic rays with the ambient matter and radiations 
during their propagation in the Milky Way has been calculated in this 
paper. If the primary ultra high energy cosmic rays are produced 
from Galactic point sources then those point sources are also emitting 
PeV gamma rays. We discuss that the detection of galactocentric PeV 
gamma rays in the future would be a signature of the presence of EeV 
cosmic accelerators in the Milky Way. 

PACS numbers: 95.85.Pw, 98.35.Jk, 98.35. Gi 

1 Introduction 

The observational data from extensive air shower experiments like Hires [1], 
AGASA [2J, Yakutsk ^ and Auger ^ have been widely used to undestand 
the origin of ultra high energy cosmic rays (UHECRs) [51 El E] • The max- 
imum energy of cosmic rays Emax from a cosmic accelerator is determined 
by the size (R), magnetic field (B) of the acceleration region, shock velocity 
and charge of the cosmic ray particle (Z), E^ax = Zef3BR fS]. It is not 
yet known what could be the maximum energy of the cosmic rays originated 
from sources inside Milky Way. Most of the cosmic ray sources in the Milky 
Way are located near its center, as a result a galactocentric anisotropy is 
expected in the UHECR spectrum. The observed UHECR spectrum does 
not show such anisotropy upto lEeV energy. The properties of cosmic rays 
in the energy range of O.lPeV and lEeV have been studied by KASCADE- 
Grande experiment [0]. It has been shown earlier that the knee like structure 
near 4 PeV in the all particle spectrum is due to the decrease in the flux of 
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light nuclei [TUl E]. The all particle energy spectrum of cosmic rays recon- 
structed from the KASCADE-Grande data using three different techniques 
and the hadronic interaction models QGSJET II, FLUKA shows concavity 
and a small discontinuity near 0.01 EeV, 0.1 EeV respectively [12]. The en- 
ergy dependence of composition of UHECRs can not be determined directly 
as their flux decreases with increasing energy. It is determined indirectly by 
observing extensive air showers in the atmosphere. The average depth of 
the shower maximum depends on the energy and mass of the primary parti- 
cle. The current observational results tell us above lEeV there is a gradual 
increase in average mass of cosmic rays upto energy 59EeV [13]. Recently 
it has been discussed fTl] that cosmic ray nuclei of energy lEeV to lOEeV 
can be originated from old gamma ray bursts or supernova explosions in our 
Galaxy and their directions may become isotropic as they spend long time 
in the micro-Gauss magnetic field of Galaxy. The inhomogeneous magnetic 
field in our Galaxy is responsible for trapping the charged particles for long 
time. The field strength has been estimated to be between 60/iG and 400/iG 
in the inner region of the Galaxy p[5j. Inside Milky Way there are super- 
nova remnants (SNRs), pulsar wind nebulae (PWN), microquasars and many 
unidentified sources of 7 rays. The TeV gamma ray spectrum from some of 
the SNRs in our Galaxy may be explained with hadronic interaction models 
and inverse Compton emission by electrons [161 113 • SNRs have long been 
identified as attractive sites for cosmic ray accelerations. In particular, near 
the Galactic center (GC) region, there is a large population of 7 ray sources. 
The high energy photon fluxes emitted by them have been observed by HESS 
dS], FERMI [^, AGILE [50], MILAGRO [21] and other gamma ray detec- 
tors. The photon flux from GC region as observed by FERMI-LAT has been 
fitted with conventional theoretical model of cosmic ray interactions with 
ambient gas medium, inverse Compton scattering of electrons and positrons 
by interstellar radiations [22]. If one considers only the protons, then pp and 
P7 interactions are the possible mechanisms of 7 ray production. The cross- 
section of pp interactions is much higher that p7 interactions ^23j. In the GC 
region the radiation and interstellar gas densities are high compared to the 
outer regions of Galaxy. The Galactic interstellar radiation field has been 
calculated by Moskalenko et al. [21]. With the average number density of 
hydrogen gas molecules n = 10cm~^ and 120cm~^ [12] and photon density 
of Galactic radiation field [23] , one finds pp process is much more important 
than p7 process in the central region of Milky Way. In these interactions 
nearly 20% of the original cosmic ray proton's energy goes to the neutral 
pion and it decays to two high energy gamma rays of equal energy. Hence, 
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we expect 100 PeV gamma rays from lEeV cosmic ray protons. 1 EeV pro- 
ton of extra-galactic origin may appear like UHECR of Galactic origin in 
the high magnetic field of GC region. If the UHECRs are protons and they 
are of Galactic origin then a galactocentric flux of secondary PeV gamma 
rays may be expected in interactions of these protons with ambient mat- 
ter. The UHECR heavy nuclei would be deflected more in Galactic magnetic 
fields compared to the protons. They would interact with ambient hydrogen 
molecules producing pions, subsequently decaying to very high energy pho- 
tons and neutrinos. Also there may be photo-disintegration of heavy nuclei 
and photo-pion production by heavy nuclei. In this paper we calculate the 
secondary PeV-EeV gamma ray flux produced by the interactions of diffuse 
UHECRs with ambient matter and radiations in the central region of Milky 
Way. The secondary PeV gamma rays/neutrinos may be anisotropic because 
of the strong magnetic field in the GC region. Moreover, if the UHECRs are 
produced by point sources in the central region of Milky Way then we expect 
to detect a fiux of galactocentric PeV gamma rays from those point sources, 
above the PeV gamma ray background calculated in this work. Thus the 
detection of PeV gamma rays would be useful to know the maximum energy 
of the cosmic rays produced by sources in the Milky Way. 



2 Pure Hadron Interactions 

The cosmic ray fiux emitted from the sources can be assumed to be a mixture 
of light and heavy nuclei. They interact with the matter molecules in the 
Milky Way during propagation producing secondary gamma rays and neu- 
trinos. The cosmic ray protons' interactions with ambient matter protons 
can be written as, pp n.^ii'^'^ + X, where is the pion multiplicity. The 
charged pions decay to high energy leptons and the neutral pions to high 
energy gamma rays. The secondary particles production in pp interactions 
has been studied earlier in [251 ESI [2Zl ESI [2S] • 

The interaction of cosmic ray nuclei with cold ambient matter leads to the 
production of neutral pions, this pion fiux has been calculated in [31]. The 
neutral pions subsequently decay to energetic gamma rays. The emissivity 
of gamma rays has been calculated in [32]. The UHE photon fiux produced 
in pp interactions and subsequent 7r° decay is 
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where = O.lEp. tesc,p and tpp ~ 10^^ /nnsec are the escape and pp inter- 
action time scales for 7r° production of cosmic ray protons respectively, uh is 
number density of ambient hydrogen molecules in the Milky Way. The cross 
section for pp interactions is included in tpp. The proton flux is contained in 
this term 

r ^ _ o /'^'^"■■""^ dNp{E^o) dE^o 

where = ApE~(^, Ap is normalisation constant and a is spectral 

index of the proton spectrum, E^^o^^in = E^ + m^o / (AE^) . The maximum 
energy of pions is the maximum energy of the cosmic ray proton/nucleon 
E-jtO ^max = E^°-^ . In pp interactions vr*^ carries on the average 20% of the 
cosmic ray proton's energy. The spectrum-weighted moments Yq, has been 
calculated from 

F„ = / x''-^f^o{x)dx. (3) 



The function ^o(x) - 8.18x^/^ [ ^^^ Jj^'^l^,^,^ ))\j^ + 

with X = ET,o/Ep. For a = 2.2, 3 we find Ya is 0.079 and 0.02 respectively. 
The gamma ray flux F^^Ap{E^) produced in pure hadronic interactions of 
cosmic ray nuclei has been calculated with this expression. 

F-f,Ap{E^) = j^^G ApiE-y)Ya, (4) 

where tesc,A is the escape time scale of cosmic ray nuclei of mass number A. 
The flux of iron nuclei can be expressed as a power law in energy similar to 
the proton flux '^^pe{Efe) _ T)^^^ hadronic interaction time scale 

of cosmic ray nuclei can be estimated with hadronic interaction cross section 
for 7r° production o"^^^^ = ?)A.QA^/^mh ^3], where A is mass number of nuclei. 
This time scale is thadr — A~^^^ /uh x lO^^sec. The expression for Gap for 
iron nuclei can be calculated similar to Gpp, except the proton flux is to be 
replaced by flux of iron nuclei per unit nucleon energy. If the proton and 
iron nuclei flux normalisation constants are equal {Ap = Ap^) then using 
EpE = 56Ep one gets 

Gap{E,) = bG-'^+'GppiE,). (5) 

The escape time scales of UHE cosmic ray protons and nuclei are not known. 
In principle one may vary the escape time scale and source spectrum of cosmic 
rays to fit the observed gamma ray flux. The ratio of the fluxes of gamma 
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rays and cosmic ray nuclei of mass number A at the same energy can be 
expressed as 

^lAp _ tesc,A '^Yg a j^i ^^^^ 
FcR,A thadr « 

The secondary nucleons produced in pure hadron processes may also in- 
teract with ambient matter and radiations to produce high energy gamma 
rays and neutrinos. 



3 Photo-Hadron Interactions 

The cosmic ray nuclei are assumed to be propagating through an isotropic 
photon background with energy e and photon density per unit energy n[e). 
The energy of a cosmic ray nucleus Etot,n = InArrin where 7^ is the Lorentz 
factor of the cosmic ray nucleons, A is mass number and m„ is the mass of 
each nucleon in the nucleus. The time scale of photo-disintegration of cosmic 
ray nuclei in the radiation field of Milky Way is tphot-d = ^ where Xa is the 
mean free path of photo-disintegration of cosmic ray nucleus |32j 

tphot-d Xa 2 Jo 7^e^ Jo 

The cross section for photo-disintegration of a nucleus of mass number A by 
a photon of energy e' in the rest frame of the nucleus is (Ja{^')- Detailed 
discussions on cross sections for all different nuclear species can be found in 
[31] • Eqn.(I7]) can be approximated as 

= " o / -^n{e). (8) 

tphot-d -'<,/27„ 

r = SMeV is the width, e'^ = 42.65A-° '^^MeV, for A > 4 is the central value 
of the GDR energy band and cto = 1.45 x lO'^'^A cm?'. The infrared radiation 
field in our Galaxy at different distances from the GC region has been calcu- 
lated by Moskalenko et al. p3]. We have approximated the infrared photon 
number density per unit energy as n{e) = eV~^cm~^ to calculate tphot-d- 
The photo-disintegration time scale and the hadronic interaction time scale 
of iron nuclei are shown in Fig.l with solid and dashed lines respectively. 
Below 100 EeV hadronic interaction is the dominant process of energy loss 
by iron nuclei. Moreover, the gamma rays produced in photo-disintegration 
of nuclei have much lower energy than the parent nuclei, hence it does not 
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Figure 1: The time scales of photo-disintegration and hadronic interactions 
of iron nuclei with ambient hydrogen molecule density uh = 10cm~^ are 
shown by solid and dashed lines respectively. 

contribute significantly to the PeV gamma ray flux. There can be photo-pion 
production by nuclei of mass number A in the radiation field of Milky Way 
Aj — > n,r7r^'° + X. The iron nuclei of energy in the range of PeV- EeV satisfy 
the threshold energy condition for Aj interactions with the x-ray background 
[35] . We find that the time scale of this interaction is much larger compared 
to the time scales of hadronic and photo-disintegration processes. We have 
calculated the time scale of photo-pion production in interactions in the 
radiation field of GC region with photon number density n(e) = -^eV'^cm'^ 
and found it is much larger compared to the pp interaction time scale with 
ambient hydrogen molecule density uh = 10cm~^. 

4 Cosmic Microwave, Galactic and Extragalac- 
tic Radio Background 

The secondary photons of energy 1 PeV satisfy the threshold energy condition 
for production of pairs with cosmic microwave background (CMB) pho- 
tons. Higher energy photons may interact with both cosmic microwave and 
radio background photons. The interaction of UHE photons {E > IPeV) 
with CMB is suppressed due to Klein Nishina suppression of interaction cross 
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section. The Galactic and extra-galactic radio backgrounds have been mea- 
sured by ARCADE 2 (Absolute Radiometer for Cosmology, Astrophysics, 
and Diffuse Emission) instrument [36] . This instrument measures the abso- 
lute temperature of the sky in search of distortions from a blackbody spec- 
trum. It operates at centimeter wavelengths between fullsky surveys at radio 
frequencies below 3 GHz and the Far Infrared Absolute Spectrophotometer 
(FIRAS) survey at frequencies above 60 GHz. The sky temperatures mea- 
sured by ARCADE 2 includes contributions from Galactic, extra-galactic 
and cosmic sources. The 2.7 K cosmic microwave background dominates 
the measured temperatures. Galactic emission is dominated by synchrotron 
emission with a smaller contribution from free-free sources. The integrated 
contribution of similar synchrotron and free-free emission in external galax- 
ies constitutes an isotropic extra-galactic radio background. Although, the 
CMB may be distinguished from Galactic or extra-galactic radio emission 
by the different frequency dependences, spectral fitting alone can not distin- 
guish Galactic emission from an extra-galactic component of similar spectral 
behavior. The method of determining the Galactic emission by using the 
north and south polar caps as convenient reference line of sights has been 
disussed in [27] . The measurement of radio emission from Galactic polar cap 
regions can be fitted with a single power law over the frequency range of 22 
MHz to 10 GHz 



with spectral index {3 = -2.55 ± 0.03 and amplitude Tq = 0.498 ± 0.028 
K at reference frequency z/q = 1 GHz. The total Galactic emission can be 
estimated with a simple plane-parallel structure of Galaxy and a cosecant fit 
in Galactic latitude angle. The ARCADE 2 data show a rise in temperature 
of 50 ± 7 mK at 3.3 GHz in addition to a CMB temperature of 2.730 ± 0.004 
K due to the extra-galactic background. The extr-agalctic radio background 
temperature has been measured in [38]. A power law fit to parametrise the 
extra-galactic background gives 



where Tcmb is the CMB baseline temperature 2.725±0.001 K, A = 1.26±0.09 
K is the power law amplitude at 1 GHz, /3 = —2.6 ± 0.04 is the power law 
index, in the frequency range of 22 MHz to 10 GHz. The secondary UHE pho- 
tons of energy 100 PeV produced in 7r° decay, satisfy the threshold condition 
for pair production with the radio photons of frequency 0.6 GHz. For the 
lower energy UHE photons energy loss due to Galactic radio background is 



TGal{l^)=To{u/uor, 



(9) 



(10) 
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not important. Radio Air-Shower Test Array (RASTA) would be an exten- 
sion of IceCube detector at the South Pole to detect air showers with surface 
array of radio antennas [39l BOj . UHE gamma rays may be detected in fu- 
ture with this detector. Photon induced showers will show a hundred times 
lower muon flux than hadron induced showers. Moreover, hadron induced air 
showers are more energetic and IceCube will have a nearly 100% efficiency of 
detecing showers of energy more than 1 PeV if the shower axis intersects the 
geometric volume of the detector. This can be used to discriminate hadron 
induced showers from photon induced showers. 

PeV photons of extra-galactic origin are unlikely to reach us due to 
their cosmological distances and higher absorption by background radiations. 
Thus the observation of PeV gamma rays would provide an exclusive window 
to study the interactions of UHECRs in the Milky Way. 

5 Optical Depths in the Galaxy 

The optical depths for very high energy photons due to background radiations 
(radio to x-rays) were calculated by Gould & Schreder [H]. The mean free 
path lenghts for interactions of UHE photons with extra-galctic radio back- 
ground down to KHz frequency was calculated by Protheroe & Biermann 
[i2] . We use the power law fits of sky temperatures due to the Galactic 
and extra-galctic radio backgrounds in the frequency range of 22 MHz to 10 
GHz given in [371 [38] to calculate the mean free path lengths for pair(e^) 
production by very high energy gamma rays. The energy density spectrum 
of Galactic and extra-galactic radio backgrounds is 

—y[BcRB{y) + BGai{y)] = ^-^{TcRB{y) + TGai{y)). (11) 

c 

The number density of radio photons per unit energy is 

dn{t) ^ 47r BcRBje) + Bcaije) 
de he e 

The mean free path of pair production has been calculated following the 
formalism discussed in ^3] . The threshold frequency of radio photons for pair 
production with 10 EeV gamma rays is nearly 6 MHz. We have assumed the 
power law fits of sky temperatures in frequency is valid upto this frequency. 
Our result for the radio background is shown in Figj2l It is comparable to 
the result given in flS]. The mean free path length for pair production due 
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Figure 2: Mean free path of pair production, dashed hne: CMB, IR from 
[24:\ : sohd hne: radio background (Galactic and extra-galactic) calculated 
using the radio background measured in [371 EH] 



to CMB is 10 Kpc for 1 PeV gamma rays |42j, which is nearly equal to the 
distance of GC from us (~ 8Kpc). For higher energy gamma rays it increases. 



6 PeV Gamma Ray Event Rates 

We find that pure hadronic interactions are more important compared to 
photohadronic processes in our study. In FigJT] the hadronic interaction 
and photo-disintegration time scales are shown for iron nuclei with dashed 
and solid lines respectively. The functions Gpp{E^), Gap{Ej) in eqn.([5]) and 
eqn.([5]) are plotted in Figj3l The upper limit on the ratio of gamma ray 
to normal cosmic ray fluxes at the same energy was set at 2.4 x 10"^ at 
310 TeV with 90% confidence limit by CASA-MIA experiment pj. The 
spectrum weighted moment in eqn.Q Ya = 0.079,0.02 for spectral index 
a = 2.2, 3 respectively. If we apply the limit from CASA-MIA experiment in 
the case of PeV gamma rays then one has to assume that the ratio of escape 
and hadronic interaction time scales of iron nuclei is ~ 0.03, 5.6 for a = 2.2, 3 
from eqn.([n]) assuming all the gamma rays are produced by hadronic interac- 
tion of UHECR iron nuclei. If we assume all the gamma rays are produced 
by UHECR protons with spectral index a = 2.2, 3 then the ratio of escape 
and pp interaction time scales of protons is 0.33 x 10~^, 1.8 x 10~^ respec- 
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tively. The high energy gamma ray event rates have been calculated assum- 
ing only iron nuclei at the source with source spectrum '^^^^^''"^ = Ap^Ep", 
Ape = IGeV^^m^'^ sec~^ sr~^ and spectral index a = 2.2. For Km^ area de- 
tector the number of events expected in one year have been plotted against 
energy, shown with solid line in FigJH We have not included the efficiency 
of the detector in this work. High energy photon absorption by CMB back- 
ground at 2 PeV energy may reduce our calculated event rate by a factor of 
0.3. At higher energies it is even less. We have multipled our calculated event 
rates by this factor to obtain a conservative estimate. As shown in Fig|2]the 
mean free path of pair production in radio background is of the order of 
Mpc, hence it is not important for these gamma rays. The observed diffuse 
UHECR spectrum has a spectral index 3 between 1 PeV and 3 EeV [15]. 
Assuming Fe nuclei and protons in 1:1 ratio we have calculated the expected 
event rates for this spectrum in Km^ — yr, shown by dotted line in Fig. 4. 
We expect a significant number of gamma ray events at PeV energy. The 
number of electron and muon neutrinos, antineutrinos produced in hadronic 
interactions is comparable to the number of gamma rays produced at the 
same energy. It is much lower than the limit set by ANTARES |46j exper- 
iment at 7.5 X 10~^{Eiy/GeV)~'^GeV~^cm~^sec~^sr~^ which corresponds to 
3 X 10^ events in Km"^ — yr. If there are point sources of very high energy 
gamma rays near the GC, PeV gamma rays may be detectable from them 
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Figure 4: High energy gamma ray event rates in Km? —yr, solid line: source 
spectrum only iron with Ap^ = lGeV~^m~'^sec~^sr~^ , a = 2.2; dashed line: 
iron and proton Ap^ — Ap — lGeV~^m~'^sec~^sr~^, a — 2.2; dotted line: 
iron and proton Ape — Ap — 1.45 x 10^ G eV ~^ m~'^ sec'^ sr"^ , a — 3. 

above the background calculated in this paper. A galactocentric very high 
energy gamma ray flux (>PeV) would be signature of EcV Galactic UHE- 
CRs. Thus the detection of PeV gamma rays would be useful to know about 
the origin of the UHECRs. 

7 Conclusions 

While propagating in the Galaxy the diffuse UHECRs interact with ambient 
matter and radiations and produce secondary particles. Some of the UHE- 
CRs may be trapped in magnetic field in the GC region and produce galac- 
tocentric secondary gamma rays. If the UHECRs are coming from Galactic 
point sources, located near the GC, then we would expect a galactocentric 
flux of gamma rays from these point sources above the background calculated 
in this work. Thus PeV gamma ray astronomy would be very useful to know 
the highest energy of the Galactic cosmic rays. 
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